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Abstract. We derive, from conformal invariance and quantum 
gravity, the multifractal spectrum /(a) of the harmonic measure (or 
electrostatic potential, or diffusion field) near any conformally invariant 
fractal in two dimensions. It gives the Hausdorff dimension of the set of 
points where the potential varies with distance r to the fractal frontier 
as r a . First examples are a random walk, i.e., a Brownian motion, 
a self-avoiding walk, or a critical percolation cluster. The generalized 
dimensions D (n) as well as the multifractal functions / (a) are derived, 
and are all identical for these three cases. The external frontiers of a 
Brownian motion and of a percolation cluster are thus identical to a 
self-avoiding walk in the scaling limit. The multifractal (MF) function 
f(a, c) of the electrostatic potential near any conformally invariant 
fractal boundary, like a critical O(N) loop or a Q -state Potts cluster, 
is solved as a function of the central charge c of the associated conformal 
field theory. The dimensions Dep of the external perimeter and Dh 
of the hull of a critical scaling curve or cluster obey the superuniversal 
duality equation (Dep — 1)(^h — 1) = \- Finally, for a conformally 
invariant scaling curve which is simple, i.e., without double points, we 
derive higher multifractal functions, like the universal function /2(a, a') 
which gives the Hausdorff dimension of the points where the potential 
jointly varies with distance r as r a on one side of the curve, and as r a 
on the other. The general case of the potential distribution between the 
branches of a star made of an arbitrary number of scaling paths is also 
treated. The results apply to critical O(N) loops, Potts clusters, and to 
the SLE K process. We present a duality between external perimeters 
of Potts clusters and O(N) loops at their critical point, as well as the 
corresponding duality in the SLE K process for kk' = 16. 

I Introduction 

Quantum field theory can be described very generally in terms of the statistics of 
Brownian paths and of their intersections 0. This equivalence is used in polymer 
theory Q| and in rigorous studies of second-order phase transitions and field the- 
ories 0] . In probability theory, non trivial properties of Brownian paths have led 
to intriguing conjectures. Mandelbrot H suggested for instance that in two dimen- 
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sions, the external frontier of a planar Brownian path has a Hausdorff dimension 
D = 4/3, identical to that of a planar self-avoiding walk, i.e., a polymer. Families 
of universal critical exponents are associated with intersection properties of sets of 
random walks PHlOf . In his famous article " Walks, walls, wetting and melting" H , 
Michael Fisher considered groups of ID "vicious walkers" which had to avoid each 
other: this amounts to considering 2D directed mutually- avoiding random walks. 
The present paper deals with the scaling properties of 2D mutually-avoiding ran- 
dom walks or fractal sets. Here the latter are conformally invariant, allowing the 
use of techniques borrowed from "2D quantum gravity" , where mutual avoidance is 
equivalent to solving a linear problem. In Q , the directness of the walks made the 
problem tractable. The general similarity in the scaling concepts used, however, 
certainly allows one to dedicate this work as a tribute to the outstanding influence 
Michael Fisher had in the general field of random scaling paths. 

The concepts of generalized dimensions and associated multifractal (Mf) mea- 
sures have been developed two decades ago |]ll|-|l4|| . Universal geometrical fractals, 
e.g., random walks, polymers, Ising or percolation models are essentially related 
to standard critical phenomena and field theory, for which conformal invariance in 
two dimensions (2D) has brought a wealth of exact results (see, e.g., ]l5|-^l[. Mul- 
tifractals and field theory must have deep connections, since the algebras of their 
respective correlation functions reveal intriguing similarities |22|. 

In classical potential theory, i.e., that of the electrostatic or diffusion field near 
random fractal boundaries such as diffusion limited aggregates (DLA), or the fractal 
structures arising in critical phenomena, the self-similarity of the latter is reflected 
in a multifractal (MF) behavior of the potential. In DLA, the potential, also called 
harmonic measure, actually determines the growth process and its scaling properties 
are intimately related to those of the of the cluster itself P3| . In statistical fractals, 
the Laplacian field is created by the random boundary, and should be derivable, in 
a probabilistic sense, from the knowledge of the latter. A first example was studied 
in Ref. [Q , where the fractal boundary, the "absorber" , was chosen to be a simple 
random walk (RW), or a self-avoiding walk (SAW), accessible to renormalization 
group methods near four dimensions. 

In two dimensions (2D), conformal field theory (CFT) has lent strong support to 
the conjecture that statistical systems at their critical point, in their scaling (con- 
tinuum) limit, produce conformally invariant (CI) fractal structures, examples of 
which are the continuum scaling limits of RW's, SAW's, critical Ising or Potts clus- 
ters (which presented a mathematical challenge, see, e.g., |25-27|). The harmonic 
measure near a such clusters possesses universal multifractal exponents, as we shall 
see. In analogy to the beautiful simplicity of the classical method of conformal 
transforms to solve 2D electrostatics of Euclidean domains, a universal solution is 
indeed possible for the planar potential near a CI fractal. 

A first exact example has been given for the whole universality class of random or 
self- avoiding walks, and percolation clusters, which all possess the same harmonic 
MF spectrum in two dimensions p8|-|3"c| (Related results can be found in ]3l]-|3^].) 
After a detailed description of this class and its link to quantum gravity, we address 
the general solution for the potential distribution near any conformal fractal in 2D 
|3~if| . The exact multifractal spectra describing the singularities of the potential, 
or, equivalently, the distribution of local "electrostatic" wedge angles along the 
fractal boundary, are given, and shown to depend only on the so-called central 
charge c, a parameter which labels the universality class of the underlying CFT. 
A further result, first obtained in |3~jf| , is the existence of a duality between the 
external frontiers of random clusters and their hulls, which applies in particular 
to Fortuin-Kasteleyn clusters in the Potts model, and to the so-called stochastic 
Lowner evolution (SLE) process (see below). 



A new feature will be the consideration of higher multifractality, which occurs in a 
natural way in the joint distribution of potential on both sides of a random CI scaling 
path, or more generally, in the distribution of potential between the branches of a 
star made of an arbitrary number of CI paths. The associated universal multifractal 
spectrum will depend on two variables, or more generally, on m variables in the case 
of an rn-arm star. We shall derive it first for Brownian motion or self-avoiding walks 
|3q ], before addressing the general case. 

Consider a two-dimensional very large "absorber" S, which can be a random 
walk, a self-avoiding walk, a percolation cluster, or, more generally, a (critical) 
scaling path or cluster. Define the harmonic measure H (w) as the probability that 
a random walker (RW) launched from infinity, first hits the outer "hull's frontier" 
or (accessible) frontier J-(S) at point w € ^F{S). One then considers a covering of 
T by balls B(w, a) of radius a, and centered at points w forming a discrete subset 
T I {a} of T . Let H(.F n B(w, a)) be the harmonic measure of the points of T in the 
ball B(w, a). We are interested in the moments of H, averaged over all realizations 
of RW's and S 

Z n = l H n (^n%a))V (1) 

\weF/{a} I 

where n can be, a priori, a real number. For very large absorbers S and hull's 
frontiers T (S) of average size R, one expects these moments to scale as 

Z n » (a/R)^ , (2) 

where the radius a serves as a microscopic cut-off, reminiscent of the lattice struc- 
ture, and where the multifractal scaling exponents r (n) encode generalized dimen- 
sions 

D(n) = ^, (3) 

71—1 



which vary in a non- linear way with n |1 1| — |X-4J . Several a priori results are known. 
D(0) is the Hausdorff dimension of the accessible frontier of the fractal. By construc- 
tion, H is a normalized probability measure, so that r(l) = 0. Makarov's theorem 



|3q| , here applied to the Holder regular curve describing the frontier |37| , gives the 
non trivial information dimension t'(1) = D (1) = 1. The multifractal formalism 
|ll] [lj] further involves characterizing subsets T a of sites of the hull's frontier T 
by a Holder exponent a, such that the H-measure of the frontier points in the ball 
B(w,a) of radius a centered at w a € T a scales as 

H(Jn% n ,a)) w (a/R) a . (4) 

The Hausdorff or "fractal dimension" / (a) of the set T a , such that 

Gard^aPBij'W, (5) 

is given by the symmetric Legendre transform of r (n) : 

a = — (n) , t (n) + / (a) = an, n = -J- (a) . (6) 
dn da 

Because of the statistical ensemble average (|l|), values of / (a) can become negative 
for some domains of a J24J. As we shall see, the associated exponents r(n) above 



can be recast as those of star copolymers made of independent RW's in a bunch, 
diffusing away from a generic point of the absorber. 

One can equivalently consider potential theory near the same fractal boundary, 
now charged. One assumes the absorber to be perfectly conducting, and introduces 
the harmonic potential H(z) at points z in the domain exterior to J 7 , with (Dirichlet) 
boundary condition H = on J 7 , and H = 1 on a large exterior circle. Then the 
local behavior of the potential 



H(z^w a ) ~r Q , r= \z - w a \ , (7) 

i 

the local behavior ((7). 



depends on the same a-exponent as the harmonic measure around point w a G J- a , 
and f(a) = dimJ- a appears as the Hausdorff dimension of boundary points inducing 




FIG. 1. Double distribution of harmonic potential H on both sides of a simple scaling 
curve (here a SAW, courtesy of Tom G. Kennedy). The local exponents on both sides of 
point w — w ai0l i are a and a'. The Hausdorff dimension of such points along the SAW is 
f 2 (a,a'). 

Generalizations to higher conformal multifractality can be defined as follows. 
When it is simple, i.e., without double points, the conformally scaling curve T can 
be reached from both sides, with a distribution of potential H + on one side, and 
H- on the other, such that 

H+(z w a , a >) ~ r a , H_{z^ w a , a ,) ~ r a ' , (8) 

when approaching a point w a , a > of subset T a ^ a i at distance r = \z — w a , a '\ (Fig. 
1). Then a double-multifractal spectrum /2(cc, a') — dim T a%a i yields the Hausdorff 
dimension of the set of points of type (a, a'). This can be generalized to the multi- 
ple scaling behavior r a \ i = 1, ...,m of the potential in the m sectors of an m-arm 
scaling star, with a multifractal spectrum f m ({cti}), to be calculated below. 



We shall use conformal tools (linked to quantum gravity) , which allow the math- 
ematical description of random walks interacting with CI fractal structures, thereby 
yielding a complete, albeit probabilistic, description of the potential. The results 
are applied directly to well-recognized universal fractals, like O(N) loops or Potts 
clusters. In particular, a subtle geometrical structure is observed in Potts clusters, 
where the external perimeter (EP), which bears the electrostatic charge, differs from 
the full cluster's hull. Its fractal dimension Dep is obtained exactly, generalizing 
the case of percolation elucidated in J38| . We obtain a duality relation: 

(D EP -l)(Z? H -l) = i > (9) 

where Z?h > £>ep is the hull dimension. Notice that the symmetric point of (^|) 
is 3/2, which is the maximum dimension of a simple (i.e., without double points) 
conformally invariant random curve in the plane. This duality, which actually gives 
clusters EP's as simple 0(n) loops at their critical point, was first obtained in J34|] . 
It predicts a corresponding duality in the SLE process (see below). 

The quantum gravity techniques used here are not yet widely known in statisti- 
cal mechanics, since they originally belonged to string or random matrix theory. 
These techniques, moreover, are not yet within the realm of rigorous mathematics. 
Contact will be made with rigorous results recently obtained in probability the- 
ory for Brownian motion and conformally invariant scaling curves [ ^lp^j39| -p[| , or 
percolation |fl2"|-|44|| , which, by completely different techniques (using in particular 
the so-called "stochastic Lowner evolution" (SLE) Q), parallel those of statistical 
mechanics and quantum gravity. In particular, our duality equation @ brings in 
the kk' = 16 duality, where the SLE K r process is the simple frontier of the SLE K , 
for k > 4. This of course hints at deep connections between probability theory 
and conformal held theory. In particular, the correspondence extensively used here, 
which exists between scaling laws in the plane, and on a random Riemann surface 
appears as fundamental. 



II Intersections of Random Walks 



Let us first define intersection exponents for random walks or Brownian motions, 
which, while simpler than the multifractal exponents considered above, in fact gen- 
erate the latter. Consider a number L of independent random walks (or Brownian 
paths) B^'jl = 1, ..,L in Z d (or R d ), starting at fixed neighboring points, and the 
probability 

P L (t) = P {u^ =1 (£« [0, t] n [0, t]) = 0} , (10) 

that the intersection of their paths up to time t is empty At large times and 

for d < 4, one expects this probability to decay as 

P L (t)^t-^, (11) 

where £l (g?) is a universal exponent depending only on L and d. Above the upper 
critical dimension d = 4, RWs almost surely do not intersect. The existence of 
exponents £l in d = 2,3 and their universality have been proven fl(i|] , and they 
can be calculated near d — 4 by renormalization theory A generalization was 
introduced Q for L walks constrained to stay in a half-plane, and starting at neigh- 
boring points on the boundary, with a non-intersection probability Pl (t) of their 
paths governed by a "surface" critical exponent £l such that 



(12) 



It was conjectured from conformal invariance arguments and numerical simula- 
tions that in 2D j| 



and for the half-plane 



KL^h[ c -l ) +2 = 1 -L{l + 2L), (141 



3 

(c) 

where h p / q denotes the Kac conformal weight 

h (c) = [(m + l)p- mqf- 1 ^ 
P:9 4m (m + 1) 

of a minimal conformal field theory of central charge c = 1 — 6/m (m + 1) , m € N* 
Jig ]. For Brownian motions c = 0, and m = 2. For L = 1, the intriguing d = 1/8 
is actually the disconnection exponent governing the probability that the origin of 
a single walk remains accessible from infinity without crossing the walk. 

To derive the conjectured intersection exponents above, the idea |2^] is to map 
the original random walk problem in the plane onto a random lattice with planar 
geometry, or, in other words, in presence of two-dimensional quantum gravity |p5| . 
The key point is that the random walk intersection exponents on the random lattice 
are related to those in the plane. Furthermore, the RW intersection problem can 
be solved in quantum gravity. Thus, the exponents C,l (Eq. (|T^)) and (Eq. ( |l4|) 
in the standard Euclidean plane are derived from this mapping to a random lattice 
or Riemann surface. 

Random surfaces, in relation to string theory Gfil, have been the subject and 
source of important developments in statistical mechanics in two-dimensions. In 
particular, the discretization of string models led to the consideration of abstract 
random lattices G, the connectivity fluctuations of which represent those of the 
metric, i.e. pure 2D quantum gravity [47j ], One can then put any 2D statistical 
model (like Ising model [Q, self-avoiding walks ]49|] ) on the random planar graph 
G, thereby obtaining a new critical behavior, corresponding to the confluence of the 
criticality of the random surface G with the critical point of the original model. The 
critical system "dressed by gravity" has a larger conformal symmetry which allowed 
Knizhnik, Polyakov, and Zamolodchikov (KPZ) [|45[|50| ] to establish the existence of 
a relation between the conformal dimensions - 1 of scaling operators in the plane 
and those in presence of gravity, A : 

A< >=A(A-7)/(l-7), (16) 

where 7 is a parameter related to the central charge of the statistical model in the 
plane: 

C =l- 6 7 2 /(l~7); (17) 

for a minimal model of the scries (|l5|), with 7 = —1/m, and Ap ^ = hp) q . 

Let us now consider as a statistical model random walks on a random graph. We 
know H that their central charge c = 0, whence m — 2, 7 = —1/2. Thus the KPZ 
relation becomes 

A<°> =U(A) = \ A(1 + 2A), (18) 



which has exactly the same analytical form as the conjecture (|l4|)! Thus, from 
the KPZ equation one infers that the planar Brownian intersection exponents Eqs. 
(13 14|) are equivalent to Brownian intersection exponents in quantum gravity: 



Ai = i(L-i), (19) 



A L = L, (20) 

Let us now sketch the derivation of the latter quantum gravity exponents pi} ] . 

Consider the set of planar random graphs G, built up with, e.g., trivalent vertices 
tied together in a random way. The topology is fixed here to be that of a sphere 
(S) or a disc (V). The partition function is defined as 

^w-Es^y-^. < 21 > 

where \ denotes the Euler characteristic x = 2 (<S) , 1 (V) ; |G| is the number of 
vertices of G, S (G) its symmetry factor. The partition sum converges for all values 
of the parameter (3 larger than some critical (3 C . At j3 — > /3+, a singularity appears 



due to the presence of infinite graphs in (21) 



Z x ((3)~(f3-(3 C ) 2 -^' M , (22) 
where 7str(x) is the string susceptibility exponent. For pure gravity as described in 



(21), the embedding dimension d — coincides with the central charge c — 0, and 
7tr(x) = 2-fx 

Now, put a set of L random walks B = {-Bjj , I = 1, L} on the random graph 
G with the special constraint that they start at the same vertex i S G, end at the 
same vertex j £ G, and have no intersections in between. We introduce the L— walk 
partition function on the random lattice p8[ : 

zm,*)= E ^r^ |G| E E zle| ' (23) 

planar G ^ ' ijEG B C) 

ij 

1=1,. ..,L 

where a fugacity z is associated with the total number \B\ = ju^-B^- 1 1 of vertices 
visited by the walks. 

We generalize this to the boundary case where G now has the topology of a disc 
and where the random walks connect two sites i and j now on the boundary dG : 

Z L VJ,(f,z)= £ e-^le-^l E E *' B| ' ( 24 ) 

disc G i,j£G B W 

ij 

l=l,...,L 

where e~@ is the fugacity associated with the boundary's length. 

The double grand canonical partition function (|2^) associated with non- 
intersecting RWs on a random lattice can be calculated exactly pg|j . The critical 
behavior of Zl (/3, z) is then obtained by taking the double scaling limit (3 — > j3 c 
(infinite random surface) and z —f z c (infinite RWs). The analysis of this singular 
behavior in terms of conformal dimensions is performed by using finite size scaling 
(FSS) H|, where one must have \B\ ~ |G|*. One obtains §§: 



Z L {p,z)~ {P-f3 c ) L ~ \G\- L . 



(25) 



Zl ( |23| ) represents a random surface with two punctures where two conformal oper- 
ators of dimension Ai are located (here two vertices of L non- intersecting RW's), 
and in a graphical way scales as 



Zl 



x\G\ 



-2A, 



(26) 



where the partition function of the two-puncture surface is the second derivative of 
Z x =i(J5) (p^). The latter two equations yield 

2Ai - 7str ( X = 2) = L, (27) 

where 7str(x = 2) = —1/2. We thus get the announced result 



A, 



(28) 



For the boundary partition function Zl ( p4|) a similar analysis can be performed 
near the triple critical point where the boundary length also diverges. The boundary 
partition function Zl corresponds to two boundary operators of dimensions A^, 
integrated over SG, on a random surface with the topology of a disc, or in graphical 
terms: 



) x \dG\ 



-2Ar 



From the exact calculation of the boundary partition function 
further equivalence to the bulk one: 



Zl/Z( 



)~Zl, 



where the equivalences hold true in terms of scaling behavior. 
^ ©, and H), and using the FSS \dG\ ~ |G| 1/2 gives 



eqs 



A, , = L. 



(29) 

one gets the 

(30) 
Comparing 

(31) 



Applying the quadratic KPZ relation ( Jig ) to A^ and A^ above yields at once the 
values in the plane R 2 , A^ = Ql (Eq. 



1), and Kf = 2( L (Eq. (0)). 
Consider now the exponents C( n ii ••; n L) — A^ ) {n{\ , as well as 2£(ni, .., ul) — 
A^ ^ {ni} , describing L mutually-avoiding bunches I = 1, .., L, each made of ni inde- 
pendent walks, i.e., mutually "transparent" f32| , with possible mutual intersections 
in a bunch. In presence of gravity each bunch will contribute its own normalized 
boundary partition function as a factor, and yield a natural generalization of (BCj) 



Z{m} 



Z{ni} 



n 



Z(n{) 



z(O) a^ov 



(32) 



to be identified with |<9G| 



-2A{n,} 



The factorization property (p2|) immediately 



implies the additivity of boundary conformal dimensions in presence of gravity 

L 

A{m,..,?u} = ^A(nO, (33) 
i=i 



where A(n) is now the boundary dimension of a single bunch of n transparent walks 
on the random surface. We know A(n) exactly since it corresponds in the standard 
plane to a trivial surface conformal dimension A(°'(n) = n. It thus suffices to invert 



(18) to get 



A(n) = [/-V) = -(V24n+ 1 - 1). (34) 



One notes the identification (|32j), on a random surface, of the bulk partition function 
with the ratio of boundary ones. In the plane, using once again the KPZ relation 



(18) for A{n/} gives the general results p8[ 

1 



ani,..,n L ) = V(x) = ^(4x 2 -l), (35) 



2C(ni, .., ? i L ) = £/(x) = ix(l + 2a;), (36) 



L L 

^ = ^{/- 1 (n i )-^ i (V24^TT-l). (37) 
l=i 1=1 

Lawler and Werner pl[ | proved by probabilistic means, using the geometrical con- 
formal invariance of Brownian motions, the existence of two functions U and V 
satisfying the structure ( |35f|3"7j ) . The quantum gravity approach here understands 
this structure in terms of linear equation ([33j), and yields the explicit functions U(x) 
and V(x) = U(\{x — |)) of (35)(|36]). The same expression of these functions has 



been finally derived in probability theory 

Let us remark that the above equations yield for ((2, l^)) describing a two-sided 
walk and L one-sided walks, all mutually non-intersecting, 

C(2, lW) = C L+i = V(L + |) = i(L + 1)(L + 2). (38) 

For L = 1, C(2, 1) = C5/2 = 1 gives correctly the escape probability of a RW from 
another RW. For L = 0, C(2, = C3/2 = 1/3 is related to the Hausdorff dimension 
of the frontier by D = 2 - 2C [§. Thus we obtain (2§] 

D = 2-2C| = | (39) 

i.e., Mandelbrot's conjecture. This conjecture has been finally proven in probability 
theory |4l| , using the analytic properties of the exponents derived from the so-called 
stochastic Lowner equation |p9|| . The quantum geometric structure explicited here 
allows several generalizations, which we now describe [ p9| . 

Ill Random Walks and Self-avoiding Walks 

We now generalize the scaling structure obtained in the preceding section to 
arbitrary sets of random or sel-avoiding walks [pi (see also J3l],[32]] ) . Consider a 
general star copolymer S in the plane R 2 (or in Zi), made of an arbitrary mixture 
of Brownian paths or RW's (set B) , and polymers or SAW's (set V) , all starting 
at neighboring points. Any pair (A, B) of such paths, A,BeB or V, can be con- 
strained in a specific way: either they avoid each other (A fl B = 0, noted A A B) , 



or they are independent, i.e., "transparent" and can cross each other (noted Ay B) 
[29|j5^| . This notation allows any nested interaction structure [p9f ; one can decide 
for instance that the branches {Pi € l of an L-star polymer, all mutually- 

avoiding, further avoid a bunch of Brownian paths {Bk G B} k=1 , all transparent 
to each other: 



5 = ( A* =1 Pt ) A (V /b= i Bk ) ■ ( 4 °) 



In 2D the order of the branches of the star copolymer does matter and is intrinsic 
to our (A, V) notation. 

To each specific star copolymer center S is attached a conformal scaling operator 
with a scaling dimension x (S) . To obtain proper scaling we consider the partition 
functions of Brownian paths and polymers having the same mean size R. When the 
star is constrained to stay in a half-plane with its core placed near the boundary, its 



partition function will scale with new boundary scaling dimension i (5) [g 19 2Q l. 

Any scaling dimension x in the bulk is twice the conformal dimension (c.dT) A^ ^ 
of the corresponding operator, while near a boundary (b.c.d.) they are identical: 

a; = 2A(°), x = hS°\ (41) 

As above, the idea is to use the representation where the RW's or SAW's are 
on a 2D random lattice, or a random Ricmann surface, i.e., in the presence of 2D 
quantum gravity ]45[ |. The general relation ( p"8| ) depends only on the central charge, 
and is valid for polymers, for which c = 0. Let us summarize the results f29| , 
expressed here in terms of the scaling dimensions in the standard plane. For a 
critical system with central charge c = 0, the two universal functions: 

U (x) = \x{l + 2x) , V (x) = ^ (4.x 2 - 1) , (42) 

with V {x) = U (2 [x — fj) i generate all the scaling exponents. The scaling expo- 
nents x(A A B), and x (A A B) , of two mutually avoiding stars A, B, with proper 
scaling exponents x (A) , x (B) , or boundary exponents x (A) , x, (B) , obey the star 



x{AAB) =2V [U- 1 (2 (A)) + U- 1 (x (B))] 

x(AAB) =U [U- 1 (x (A)) + If- 1 (x {B))] , (43) 
where U~ x (x) is the inverse function of U 

U^ 1 (x) = - (V24x + 1 - 1) . (44) 

On a random surface, U~ x (x) is the boundary dimension corresponding to the 
value x in R x R + , and the sum of U^ 1 functions in Eq. ( |4^ ) represents linearly 
the juxtaposition A AB of two sets of random paths near their random frontier, i.e., 
the product of two "boundary operators" on the random surface. The latter sum 
is mapped by the functions U, V, into the scaling dimensions in R 2 psfl . 

The rules (p!3|), which mix bulk and boundary exponents, come from simple fac- 
torization properties on a random Riemann surface, i.e., in quantum gravity [2^j29[| , 
(and are also recurrence relations in R 2 between conformal Riemann maps of the 
successive mutually-avoiding paths onto the line R [^H). 

If, on the contrary, A and B are independent and can overlap, then by trivial 
factorization of probabilities their dimensions are additive |29] 



x(AVB) = x (A) + x (B) , 
x{AV B)=x {A) +x(B). 



(45) 



It is clear at this stage that the set of equations above is complete. It allows for 
the calculation of any conformal dimensions associated with a star structure S of 
the most general type, as in (fiO|), involving (A,V) operations separated by nested 
parentheses p9| . 

Brownian-polymer exponents: The single extremity scaling dimensions are for a 



RW or a SAW near a Dirichlet boundary in R 2 [EO 55 



x B (1) = A<°> (1) = 1, x P (1) = A^ (1) = |, (46) 

or on G, A B (1) = U~ x (1) = 1, A P (1) = U^ 1 (§) = f . Because of the star algebra 
described above these are the only numerical seeds, i.e., generators, we need. 

Stars can include bunches of n copies of transparent RW's or m transparent 
SAW's. Their b.c.d.'s in R 2 are respectively, by using fl45| ) and (f46|), A^ (n) = n 
and Ap^ (m) = |m, from which the inverse mapping to the random surface yields 
As (n) = U^ 1 (n) and Ap (m) = U^ 1 ■ The star made of L bunches I € 

{1, L}, each of them made of ng transparent RW's and of mg transparent SAW's, 
and the L bunches being mutually-avoiding, has planar scaling dimensions 

£(°) {n e ,m e } = U (A) , A<°> {n,,m,} = V (A) , 

A {n^, m^} = C/ _1 (nc + |m^) . 

This encompasses all previously known exponents for RW's and SAW's [p pS|j20[ |. 
We in particular arrive at the striking scaling equivalence: a self-avoiding walk is 
exactly equivalent to 5/8 of a Brownian motion. Similar results have been obtained 
in probability theory, based on the general structure of "completely conformally 
invariant processes" , which correspond exactly to c = central charge conformal 
field theories p2| , pp| . The rigorous construction of the scaling limit of SAW still 
eludes a rigorous approach, eventhough it is predicted that it corresponds to the 
"stochastic Lowner evolution" SLE K with k = 8/3, equivalent to a Coulomb gas 
with g = 4/k = 3/2 (see below section XI). 

IV Conformal Multifractality and the Harmonic Measure 

The harmonic measure, i.e., the diffusion or electrostatic field near an equipo- 
tential fractal boundary |5q] , or, equivalently, the electric charge appearing on the 
frontier of a perfectly conducting fractal, possesses a self-similarity property, which 
is reflected in a multifractal (Mf) behavior. Cates and Witten [M considered the 
case of the Laplacian diffusion field near a simple random walk, or near a self- 
avoiding walk. The associated exponents can be recast as those of star copolymers 
made of a bunch of independent RW's diffusing away from a generic point of the 
absorber. The exact solution to this problem in two dimensions is as follows p9| . 
From a mathematical point of view, it could also be derived from the results of refs 

HHSIl taken alt °g etner - 

The two-dimensional "absorber" S can be a random walk, or a self-avoiding walk. 
The harmonic measure H (w) is the probability that another random walker (RW) 
launched from infinity, first hits the outer "hull's frontier" or (accessible) frontier 
!F(S) at point w E ^F(S). A covering of T by balls B(w, a) of radius a is centered 
at points w € JF/{a} forming a discrete subset JF/{a} of T . Let H(.F D B(w, a)) be 



the harmonic measure of the intersection set between T and the ball B(w, a). The 
moments of H, averaged over all realizations of RW's and iS are defined as 




(TnB(w,a)) ), (47) 



where n can be, a priori, a real number. In the limit of large absorbers S and 
frontiers !F (S) of average size R, or small covering radius a, i.e, a/R — > 0, these 
moments scale as 

Z n « (a/P) r(n) , (48) 

where the multifractal scaling exponents r (n) encode generalized dimensions D (n) , 
r (n) = (n — 1) D (n) , which vary in a non-linear way with n [[ll|-[l4j. 

As explained in the introduction, the harmonic multifractal spectrum /(a) (Eqs 
(||J|,^)) is derived as a Legendre transform of the r(n) function. (The existence of 
the harmonic multifractal spectrum /(a) for a Brownian path has been rigorously 
established in @.) 

By the very definition of the H-measure, n independent RW's diffusing away 
from the absorber give a geometric representation of the n th moment H n , for n 
integer, and convexity arguments give the whole continuation to real values. When 
the absorber is a RW or a SAW of size R, the site average of its moments W 1 
is represented by a copolymer star partition function Zfi(S^n), where we have 
introduced the short-hand notation S^n = S A (VB) describing the co polym er star 
made by the absorber S hit by the bunch (VB) n at the apex only |gj,^9[. More 
precisely one has 

Z n a R 2 Z R (S A n) (49) 

where the absorber S is either the two-RW star B V B or the two-SAW star PAP, 
made of two non-intersecting SAWs. Owing to Eq.(p8|), we get the scaling relation 

T{n) = x (5 A n) - 2. (50) 

Our formalism (|43| ) immediately gives the scaling dimensions 

x{S A n) = 2V^A(S) + U- 1 (n)^j, (51) 

where A (S) is as above the quantum gravity boundary dimension of the absorber 
S alone. For a RW absorber, we have A(PVP) = U^ 1 (2) = |, while for a 
SAW A (P A P) = 2A P (1) = 21J- 1 (§)=§. The coincidence of these two values 
tells us that in 2D the harmonic multifractal spectra f (a) of a random walk or a 
self-avoiding walk are identical. The calculation gives 123] 

1 



r(n) = -(n-l) + — (a/24^+1-5), (52) 



dr , . 1 5 1 ,_„. 
» = o + o-7S5=tt=' ( 53 ) 
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D{ ^ = \ + ^xJr^ I r ' ne[-S'+ 00 )' ( 54 ) 
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FIG. 2. Harmonic multifractal dimensions r(n) and spectrum /(a) of a two-dimensional 
RW or SAW. 



The corresponding universal curves are shown in Fig.l: r (n) is half a parabola, 
and / (a) a hyperbola. £> (1) = r' (1) = 1 is Makarov's theorem. The singularity 
at a = i in the multifractal functions /(a) corresponds to points on the fractal 
boundary T where the latter has the local geometry of a needle. The mathematical 
version of this statement is given by Beurling's theorem [Q, which states that at 
distance e from the boundary, the harmonic measure is bounded above by 

H(z : mf w€ r\z-w\ < e)) < Ce 1/2 , (56) 

where C is a constant. This insures that the spectrum of multifractal Holder ex- 
ponents a is bounded below by b. The right branch of the/ (a) curve has a linear 
asymptote 

lim ^/(«) = -^- (57) 

Its linear shape is quite reminiscent of that of the multifractal function of the growth 
probability as in the case of a 2D DL A cluster JsTJ] . The domain of large values of a 
corresponds to the lowest part n — > n* = — ^ of the spectrum of dimensions, which 
is dominated by almost inaccessible sites, and the existence of a linear asymptote to 
the multifractal function / implies a peculiar behavior for the number of those sites 
in a lattice setting. Indeed define Af (H) as the number of sites having a probability 
H to be hit: 

J\f (H) = Card {w G T : H(to) = } . (58) 



Using the Mf formalism to change from variable H to a (at fixed value of a/R), 
shows that AT (H) obeys, for H — ► 0, a power law behavior 

N(H)\ H ^o~H- T * (59) 

with an exponent 

t* = 1+ lim -f (a) = l + n*. (60) 

Thus we predict 



23 
24' 



(61) 



One remarks that — r (0) = sup,-,, / (a) = f (3) = | is the Hausdorff dimension of 
the Brownian frontier or of a SAW. Thus Mandelbrot's classical conjecture identify- 
ing the latter two is derived and generalized to the whole f (a) harmonic spectrum. 

- An Invariance Property of f(a) 

The expression of /(a) simplifies if one considers the combination: 

f( a )-a= — [l- - ( 2a -1 + —?—)] . (62) 
J y ' 24 2 V 2a- 1 J V ' 



Thus the multifractal function possesses the invariance symmetry 60 



f(a)-a = f(a')-a', (63) 
for a and a' satisfying the duality relation: 

(2a- l)(2a'- 1) = 1, (64) 

or, equivalently 

a~ 1 +a'~ 1 =2. (65) 

When associating a wedge angle 9 = ir/a to each local singularity exponent a, one 
recovers the complementary rule for angles in the plane |30[ 

9 + 6' = - + —, = 2tt. (66) 
a a 



It is interesting to note that, owing to the explicit forms ( |53j ) of a and ( |54| ) of D(n), 
the condition (|65|) reads also after a little algebra 

D(n) + D(ri) = 2. (67) 

This basic symmetry ( |63| ) reflects that of the cluster boundary itself under the ex- 
change of interior and exterior domains (ref. ]60[|). 

- Higher multifractality of Brownian motion and self-avoiding walk 

It is interesting to note that one can also define higher multifractal spectra as 
those depending on several a variables f^l . A first exemple is given by the double 
moments of the harmonic measure on both sides of a random fractal, taken here 



as a Brownian motion or a self-avoiding walk. (The general case will be further 
described below in section VI.) The fractal boundary has to be reached from both 
sides, so it must be a simple curve without double points, which is naturally the 
case of a SAW. For a Brownian motion, one can consider the subset of the pinching 
or cut points, of Hausdorff dimension D = 2 — 2(2 = 3/4, where the path splits into 
two non-intersecting parts. Locally the Brownian path then is accessible from both 
directions. 
Let us define: 



Z n ,n' = ( £ [H+MHH- 
\w£F/{a} 



(68) 



where H + (u>) = H + (f n6(i»,o)) and H_(w) = H (J 7 C\ B(w, a)) are respectively 
the harmonic measures on "left" or "right" sides of the random fractal. These mo- 
ments have a multifractal scaling behavior 



Z n » (a/R) 



T2(n,n') 



(69) 



where the exponents r 2 (n, n') now depend on two moment orders n and n' . The 
generalization of the Legendre transform Eq. (||) reads 

a = — (n,n), a = — (n,n), 
on an' 

/2(a,a) = an + an — T 2 {n, n'), 

n = — (a, a ) , n = — (a, a ) . 

We find the r exponents from the star algebra ([43|): 

T 2 (n, n') = 2V (a' + U' 1 (n) + l/' 1 (n')) - 2, 



(70) 



(71) 



where a' corresponds to the quantum gravity scaling dimension of the fractal set, 
i.e., the simple curve or the pinching point set, where the harmonic measure is 
evaluated on both sides. For a Brownian motion, pinched into two parts separated 
by the two sets of auxiliary Brownian motions, representing the moments of the 
harmonic mesures, we so have: 

a' = A (B A B) = 2 x A B (1) = 2C/" 1 (1) = 2. (72) 

For a self-avoiding walk made of two mutually-avoiding one-sided arms, we have 



a' = A(PAP) = 2x Ap(l) = 2U- 1 ( I) = |« 



(73) 



After performing the double Legendre transform and some calculations, we find 

-1 



1 1 2 

h {a, a) = 2 + — - -a 
- 54 (" + « ), 



1 - 



1 / 1 



1 



(74) 



a = 2- 



1 

' V24n + 1 



1 + V24n' + 1 



(75) 



and a similar symmetric equation for a' . Here a" has the shifted value: 



a" = a' + 7 = a' - - (76) 

= ^ (RW), or a" = 1 (SAW). (77) 

This doubly multifractal spectrum possesses the requested properties, like 
sup Q //(a, a') — f{a), where f(a) is (|5|) above. 

This can be generalized to a star configuration made of m random walks or m 
self-avoiding walks, where one looks at the simultaneous behavior of the potential 
in each sector between the arms of the star (see VI below for a more precise 
description in the general case). The poly-multifractal results read for Brownian 
motions or self-avoiding polymers: 

1 1 / 1 m \ 1 

fm ({«i=l,..., m }) = 2 + - - -o" 2 l--^^ 1 



12 3 \ 2 

\ 2=1 



Tfl 
8=1 



(78) 



with 



2 ;ra(°" + ig^ TT )' (79) 



and where 



1 3 

a = m, a" = a' — -m = -m, (80) 

for m random walks in a star configuration, and 

3 11 

a' = — m, a" = a' m = —m, (81) 

4 4 2 

for m self-avoiding walks in a star configuration. The two-sided case above (^) is 
recovered for m = 2. The domain of definition of the poly-multifractal function / 
is given by 

i=l 

as verified by Eq. ([79]). 
V Percolation Clusters 



Consider now a two-dimensional very large incipient cluster C, at the percolation 
threshold p c . Define H (w) as the probability that a random walker (RW) launched 
from infinity, first hits the outer (accessible) percolation hull's frontier J-(C) at point 
w € J~{C). The moments of H are averaged over all realizations of RWs and C, 
as in Eq.(^7|) above. For very large clusters C and frontiers T (C) of average size 
i?, one expects again these moments to scale as in Eq. p^): Z n w {a/R) T<yn \ 
These exponents r(n) have been obtained recently |]3(i[| , using an exact result on the 
external boundary of a percolation cluster |3£| . 



We consider site percolation on the 2D triangular lattice. Figure 2 depicts n 
independent random walks, in a bunch, first hitting the external frontier of a per- 
colation cluster at a site w = (•) . This site, to belong to the accessible part of the 
hull, must remain, in the continuous scaling limit, the source of at least three non- 
intersecting crossing paths, noted £3, reaching to a (large) distance R (Fig. 3) pa ]. 
(Notice that the definition of the standard hull requires only a pair of dual lines) . 
The n independent RW's, or Brownian paths B in the scaling limit, in a bunch 
noted (\/B) n , avoid the set S3 = (AT 3 ) of three non-intersecting connected paths 
in the percolation system, and this system is governed by a new critical exponent 
x (S3 A n) depending on n. 

0-0 0-0 




FIG. 3. An accessible site (•) on the external perimeter for site percolation on the trian- 
gular lattice. It is denned by the existence, in the scaling limit, of three non-intersecting, 
and connected paths 1S3 (dotted lines), one on the incipient cluster, the other two on the 
dual empty sites. The entrances of fjords close in the scaling limit. Point (•) is first 
reached by three independent RW's (red, green, blue), contributing to H 3 («). The hull of 
the incipient cluster (golden line) avoids the outer frontier of the RW's (thick blue line). 



In terms of these definitions, the harmonic measure moments simply scale with 
an exponent E^] 



r (n) = x (S3 A n) 



(83) 



For percolation, two values of half-plane crossing exponents xi are known by ele- 
mentary means: x,i = \,x$ = 2 |2q . We fuse the two objects £3 and (\/B) n into a 
new star S3 A (yB) n , and use (|43|) to obtain 



x (S3 A n) = 2V (U- 1 (z 3 ) + U- 1 (n)) 
Specifying U' 1 (x 3 ) = § finally gives from ©(El) 



(84) 



x (S3 A n) = 2 + - (n - 1) + — (V24n + 1 - 5) 



24 

From this r (n) ( |S3"|) is found to be identical to 
follows as: 



§ for RW's and SAWs; D (n) 



valid for all values of moment order n, n > 
again exactly as in Eq. (|55 



-j. The Legendre transform reads 



/(a) = 48 



1 



2a - 1 



a 
24 



(86) 



Only in the case of percolation has the harmonic measure been systematically 
studied numerically, by Meakin et al. J6lj . We give in Figure 3 the exact curve 
D (n) ( |85| ) 30 together with the numerical results for n e {2, ...,9} [|Sl| , showing 
fairly good agreement. 




FIG. 4. Universal generalized dimensions D(n) as a function of n, corresponding to the 
harmonic measure near a percolation cluster, or to self-avoiding or random walks, and 
comparison with the numerical data obtained by Meakin et al. (1988) for percolation. 

The average number Af(H) ( |59| ) has been also determined numerically for percola- 
tion clusters in |52l, and for c = 0, our prediction (^lj) r* = || = 0.95833... compares 
very well with the result t* = 0.951 ± 0.030, obtained for 10~ 5 < H < 10" 4 . 

The dimension of the support of the measure D(0) = | ^ D K , where £> H = \ 
is the Hausdorff dimension of the standard hull, i.e., the outer boundary of critical 
percolating clusters p3{ . The value D(0) = | corresponds to dimension of the ac- 
cessible external perimeter. A direct derivation of its exact value is given in J38| . 
The complement of the accessible perimeter in the hull is made of deep fjords, which 
do close in the scaling limit and are not probed by the harmonic measure. This 
is in agreement with the instability phenomenon observed on a lattice for the hull 
dimension p4| . A striking fact is the complete identity of the multifractal spectrum 
for percolation to the corresponding results, Eqs.( p2|f35| ), both for random walks and 
self-avoiding walks. Seen from outside, these three fractal simple curves are not dis- 
tinguished by the harmonic measure, in fact they are the same, and one of the 
main conclusions of this study is that the external frontiers of a planar Brownian 
motion, or of a critical percolation cluster are identical to a critical self-avoiding 
walk, with a Hausdorff dimension D = | . As we have seen, this fact is linked to the 
presence of a single universal conformal field theory (with a vanishing central charge 
c = 0), and to the underlying presence of quantum gravity, which structures the 
associated conformal dimensions. Note that in a recent work, Smirnov ]42] | proved 
that critical site percolation on the triangular lattice has a conformally invariant 
scaling limit, and that the discrete cluster interfaces (hulls) converge to the same 



stochastic Lowner evolution process as the one involved for Brownian paths, open- 
ing the way to a rigorous derivation of percolation exponents fj^Jfil , previously 
derived in physics |l7|]55|,|63|,|38[ . 



VI General Conformally Scaling Curves and Higher Multi- 
fractality 

In the next sections, we present the main description of multifractal functions 
results in a universal way. For multiple simple curves, we also define the higher 
multifractal spectra, depending on an arbitrary number of a- variables. We then 
proceed with their derivation from conformal field theory and quantum gravity. 
The geometrical findings are described in details, including well-known cases (Ising 
clusters, Q = 4 Potts model). Finally, some geometric duality properties for the 
external boundaries in O(N) and Potts models are explained. We also make explicit 
the relation between a conformally invariant scaling curve with CFT central charge 
c ]34j , and the stochastic Lowner process SLE K p9|- ^l| , [74|j7l| ] . 

Consider a single (conformally invariant) critical random cluster, generically 
called C. Let H (z) be the potential at exterior point z € C, with Dirichlct boundary 
conditions H (w £ dC) = on the outer (simply connected) boundary dC of C, (or 
frontier T = dC) , and H(w) — 1 on a circle "at oo", i.e., of a large radius scaling 
like the average size R of C. From the well-known theorem due to Kakutani p5[ , 
H (z) is identical to the harmonic measure, i.e, the probability that a random walker 
(more precisely, a Brownian motion) launched from z, escapes to oo without having 
hit C. The multifractal formalism [pl]-^4[ characterizes subsets dC a of boundary 
sites by a Holder exponent a, and a Hausdorff dimension / (a) = dim (dC a ), such 
that their potential locally scales as 

H(z^wedC a )& {\z-w\/R) a , (87) 

in the scaling limit a«r= \z — w\ <C R, with a the underlying lattice constant. In 
2D the complex potential ip(z) (such that the electrostatic potential H(z) = ^ftip(z) 
and field |E(z)| = reads for a wedge of angle 9, centered at w: 

(p(z) = {z-w) 7l/e . (88) 

By Eq. j87|) a Holder exponent a thus defines a local equivalent "electrostatic" 
angle 9 = ir/a, and the MF dimension f(ff) of the boundary subset with such 9 is 

f(6) = f{a = n/6). (89) 
Of special interest are the moments of H, averaged over all realizations of C 



Z n = ( V H n (z) , (90) 




where dC{r) is shifted a distance r outwards from dC, and where n can be a real 
number. In the scaling limit, one expects these moments to scale as 

Z n « (r/i?) r(n) , (91) 

where the multifractal scaling exponents r (n) encode generalized dimensions, 
D (n) — t (n) I (n — 1), which vary in a non-linear way with n P1"|l4]|; they obey 
the symmetric Legendre transform r (n) + / (a) = an, with n — f (a) , a = r' (n). 



From Gauss' theorem E4] r(l) = 0. As said above, because of the ensemble average 



fl90D , values of / (a) can become negative for some domains of a |22|,|2J . 

Now, we consider the specific case where the fractal set C is a (conformally invari- 
ant) simple scaling curve, that is, it does not contain double points. The frontier 
dC is thus identical with the set itself: 

dC = C. (92) 

Each point of the curve can then be reached from infinity, and we can address the 
more refined question of the simultaneous behavior of the potential on both sides 
of the curve. Specifically, the potential H scales as 

H+ (z -> w+ e dC a , a >) w \z~w\ a , (93) 

when approaching w on one side of the scaling curve, together with the scaling 

H- w~ e dC a , a ,) k, \ z - w \ a \ (94) 

on the other side. We can then generalize the multifractal formalism to characterize 
subsets C a>a > of boundary sites w by two Holder exponents a, a' such that the 
potential near w locally scales on the two sides of C as in Eqs. ( |93| ) and (pif). 
This subset is characterized by a Hausdorff dimension f 2 (a, a') = dim (C a , a >). The 
standard multifractal spectrum f(a) is then recovered as the supremum: 

f(a) = sup a '/ 2 (a, a') . (95) 

As above, one can also define two equivalent "electrostatic" angles from the Holder 
exponents a, a', as 9 = ir/a,8' = n/a' and the MF dimension f 2(0,0') of the 
boundary subset with such 0, 8' is then 

h(0,6') = f2(a = TT/e ) a' = Tr/e'), (96) 

Define the harmonic measure H (w) as the probability that a random walker (RW) 
launched from infinity, first hits the frontier C at point w E C. A covering of C by 
balls B(w,r) of radius r is centered at points w £ C/{r} forming a discrete subset 
C/{r} of C. Let H(C n B(w,r)) be the harmonic measure of the intersection of C 
and the ball B(w, r). The double multifractal spectrum will be computed from the 
double moments of the harmonic measure on both sides of the random fractal curve. 
Let us define: 



H + H]"[H_H]" , (97) 



where H + (w) = H + (Cn8(w,r)) and H_(iy) = H_ (CC\B(w, r)) are respectively 
the harmonic measures on the "left" or "right" sides of the random fractal. These 
double moments have a multifractal scaling behavior 

Zn t n> « (r/RY 2 ( " ,n,) , (98) 

where the exponent T2 (n, n') now depends on two moment orders n,n'. The gen- 
eralization of the usual Legendre transform of multifractal formalism Eq. (^|) reads 

a = _ (n)n)) a =-^j(n,n), 
/2 (a, a') — an + an — t 2 (n, n) , 

n=—(a,a), n= — (a,a). (99) 




From definition (j97|) and Eq. (j98|), we recover the one-sided multifractal spectrum 
as 



r(n) =T 2 (n,n' = 0). (100) 
Putting the value n' = in the Legendre transform Eq. we obtain the identity 



(95), as it must. 




More generally, one can consider a star configuration S m of a number to, to > 2, 
of similar simple scaling paths, all originating at the same vertex w. The higher 
moments -Z„ liJl2j ... inm can then be defined as 

2 ni ,n 2 ,..., nm = (y, PiMP [H 2 MP • • • [H m Mp> ) , (101) 
where 

H,H = H, (Cn%,r)) 

is the harmonic measure (or, equivalently, local potential at distance r) in the zth 
sector of radius located between paths i and i + 1, with i = 1, • • • ,to, and by 
periodicity m + 1 = 1 . These higher moments have a multifractal scaling behavior 

Z nun2 ,..., nm « (r/ii)^'™*'-'™".) , (102) 

where the exponent r m (ni,n2, ...,n m ) now depends on the set of moment orders 
ni,n2, ...,n m . The generalization of the usual Legendre transform of multifractal 
formalism Eq. (||) now involves a higher multifractal function f m (at, ct2, • ■ ■ , cn m ), 
depending on to local exponents a,: 

m 

/m ({«i}) = X! aini ~ Tm ({ ni }) ' 
t=l 

^=^(K})- ( 103 ) 

At this point, a caveat is in order. The reader could wonder about the meaning 



of the sum over point w in (101), if there is only one such to- vertex in a star! This 
formal notation is kept in continuity with the to = 2 case, and annonces the fact 
that exponents r m (ni, Ti2, n m ) are calculated with inclusion of the Hausdorff 
dimension D mi associated with the star center, which reads in these notations 
D m = —T m (0, 0, 0) = sup{ a . yf m ({ctfj}), and becomes negative for to high enough 
(see section VIII below). One can define shifted exponents f m = r m + D m — 
T m [ni,n2, ...,n m ) — T m (0,0, ...,0), which correspond to a different normalization, 
and describe the scaling of local averages 

([HiHr [H 2 HP • • • [H m H]" m ) « {r / K) ^(m.n„...,n m ) _ (1Q4) 



By Legendre transform ( |103[ ) these exponents give the subtracted spectrum 
/m ({ a «}) ~ su P{ai}/m ({ a i}) directly. The latter has an immediate physical mean- 
ing: the probability P({cti}) to find a set of local singularity exponents {cki} in the 
to sectors of an m-arm star scales as: 



P m ({aJ) oc fl/»««i})/JJ»up/m . (105) 



From definition (101) and Eq. ( |102| ), we get the lower (m — l)-multifractal spec- 
trum as 

r ir~ 11 (ni,n2, ■ ■ ■ , n m -i) = r m (ni,n 2 , • • • , n m -i,n m = 0) . (106) 

In these exponents, the subscript m stays unchanged since it counts the number of 
arms of the star, while the potential is evaluated only at m — 1 sectors among the 
m possible. More generally, one can define exponents 



(ni,n 2 , ■ ■ ■ ,n p ) = r m (n\,n 2 , ■ ■ ■ ,n p ; n p _i = 0, • • • ,n m = 0) , 



where p takes any value in 1 < p < m. Note that according to the commutativity 
of the star algebra for exponents between mutually-avoiding paths (see Eq. ( fi"3| ) and 
below), the result does not depend on the choice of the p sectors among m. Putting 
the value n m = in the Legendre transform Eq. (|103|), we obtain the identity: 



/J™ _11 (ai, • • • ,a m _i) = sup Qm / m (ai,a 2 , • • ■ ,a m ) . (107) 

Note that the usual f(a) spectrum is in these notations f^\a). 

VII Conformal Invariance and Quantum Gravity 

Let us now follow the main lines of the derivation of exponents r (n) , hence 
/(a), or, more generally, r m (ni,n 2 , n m ) and f m (o>i, a 2 , ■ ■ • , a m ) by generalized 
conformal invariance. By definition of the H-measure, n independent RW's, or 
Brownian paths B in the scaling limit, starting at the same point a distance r 
away from the cluster's hull's frontier dC, and diffusing without hitting dC, give a 
geometric representation of the n th moment, H n , in Eq.(^) for n integer. Convexity 
yields analytic continuation for arbitrary n's. Let us introduce the notation AAB for 
two random sets conditioned to traverse, without mutual intersection, the annulus 
T> (r, R) from the inner boundary circle of radius r to the outer one at distance 
R, and A\J B for two independent, thus possibly intersecting, sets (29|. With 
this notation, the "probability" (actually the associated grand canonical partition 
function) that the Brownian paths and cluster are in a configuration dC A (V0) = 
dC A n, is expected to scale for R/r — > oo as 

V R (dCAn) w {r/R) x(n \ (108) 

where the scaling exponent x ( n) depends on n. In terms of definition (| 1 08| ) , the 
harmonic measure moments (|9C]) simply scale as Z n w R 2 Vr (dC A n) p^ , ^9| , which, 
combined with Eq. (^), leads to 

t (n) = x (n) - 2. (109) 

To calculate these exponents, we use the fundamental mapping of the conformal 
field theory in the plane R 2 , describing a critical statistical system, to the CFT on 
a fluctuating abstract random Riemann surface, i.e., in presence of quantum gravity 
|4|j5^,^]. Two universal functions U, and V, which now depend on the central 
charge c of the CFT, describe this map: 

U(x)=x^, V(x) = j^^, (110) 
1—7 41—7 



with 



V(x) = U{\{x + 1 )) 



(111) 



The parameter 7 is the string susceptibility exponent of the random 2D surface (of 
genus zero), bearing the CFT of central charge c fl45fl; 7 is the solution of 



1 -67 2 (1 -7)^,7 < 0. 



(112) 



For two arbitrary random sets A, B, with boundary scalin g exp onents in the half- 
plane x ( A) , x (B) , the scaling exponent x(AAB), as in (108), has the universal 
structure p9,Bq,p3 



x (A A B) =2V [U- 1 (x (A)) + U' 1 (x (B))] , 
where U~ x (x) is the positive inverse function of U 

U- 1 (x) = X - ( v /4(l- 7 )x + 7 2 
Note that one has the shift relation 

U- l {x) = \v- l {x)+ 1 ~~i 

where 







V- 1 (x) = ^4(l-7)^ + 7 2 



(113) 



(114) 



(115) 



(116) 



XJ~ X (x) is, on the random Riemann surface, the boundary scaling dimension corre- 
sponding to x in the half-plane R x R + , and the sum of t/ _1 functions in Eq. (113) 
is a linear representation of the product of two "boundary operators" on the random 
surface, as the condition A A B for two mutually-avoiding sets is purely topological 
there. The sum is mapped back by the function V into the scaling dimensions in 

R 2 M. 



For the harmonic exponents x(n) = x(dCAn) in (108), we use (113). The 
external boundary exponent x (dC) obeys 



U- 1 {x) = l- 1 , 
which we derive either directly, or from Makarov's theorem: 

dx 



r'(n = 1) 



dn 



(n = 1) = 1. 



(117) 



(118) 



The bunch of n independent Brownian paths have simply x((\/B) n ) = n, since 
x(B) = l @. Thus we obtain 



x{n) = 2V (l - 7 + If- 1 (n)) 



This finally gives from (|TTC|) (JlT4|) r(n) = x(n) - 2g: 

r (n) = i(n - 1) + ^[v/4(l- 7 )« + 7 2 - (2 - 7)] 



(119) 



(120) 



Similar exponents associated with moments later appeared in the context of the 
SLE process (see II in EJ, |fri| ) ; see also for Laplacian random walks. 



The Legendre transform is easily performed to yield: 



dr 1 
ft \ !( 2 -^ 2 

/(") = o- 



2-7 



2 ^4(1-7)71 + 72 



3- 



1 



2a- 1 



1 7 2 



1-7 
a £ (i, +00) . 

Using the identities in terms of central charge c: 

l(2 - 7 ) 2 25 -c 



41 



-a, 



(121) 
(122) 



4 I-7 
1 7 2 



41 



7 



24 
1 - c 
24 ' 



(123) 



we find 



r (n) = — (n — 1) + 



25 



24 



'24n + 1 - 
25- c 



(124) 



/(a 



n S 



25 -c 
48 

a G (i, +00 



1 



c 



-00 



24 ' 

1 1 

2 + 2 V 24// 
1 



25 -c 



1 



3- 



2a- 1 



- c 
c 



24 



(125) 



This formalism immediately allows generalizations. For instance, in place of n ran- 
dom walks, one can consider a set of n independent self-avoiding walks V, which 
avoid the cluster fractal boundary, except for their common an chori ng point. The 
associated multifractal exponents x (dC A (VP)") are given by (119), with the ar- 



gument n in U (n) simply replaced by i((VP) ) = nx (V) 



29. These 



exponents govern the universal multifractal behavior of the moments of the prob- 
ability that a SAW escapes from C. One then gets a spectrum / such that 
/ (a = x (V) 7r/0) = f (a — w/0) = f(0), thus unveiling the same invariant un- 



derlying wedge distribution as the harmonic measure, (see also 1 33 
VIII Higher Multifractal Spectra 



In analogy to Eqs. ( 109 ), (119), the exponent T2(n, n') is associated with a scaling 
dimension ^(n, n') by 



T 2 {n,n') = x 2 (n,n') - 2 

x 2 (n, n) = 2V [l - 7 + U' 1 (n) + U^ 1 (n) 



(126) 



The calculation of the double Legendre transform Eq. Q9Sj) is as follows. We start 
with the notation for the total quantum gravity scaling dimension: 



5= l--f + U- 1 (n) + U- 1 (n') 



(127) 



This gives explicitly: 



6=1 + -\/4(l - 7)71 + 7 2 + -a/4(1 - 7)"' + 7 2 



(128) 



Then, we have 



and since 



we finally get 



dx2 1 >\ oi/'m 9,5 



21-7' 



(129) 



V4(l -7)« + 7 2 ' 
A useful consequence is the identity 



v/4(l - 7)71' + 7 2 



(130) 



1 / 1 1 



2\a a' 



Equation ( 130 ) can be inverted into 



4(1-7) V« 2 



a 



(131) 



(132) 



(133) 



Reordering the S terms with use of (110), and recalling identity ( |13l| ) for 5, finally 
gives after some calculations the explicit formulae 



where use was made of ( |110[ ) for V. This allows the simple expression of f 2 
h {a, a') =2- V(6) + aV (^j + a'V (J^j . 



f 2 (a, a') = 



25 -c 



1 



12 2(1-7) 
1 - c 
~24~ 



1/1 1 
1-r - + — 



2 V a a' 



(a + a') , 



(134) 



^4(1-7)7^ + 72 



+ o (^4(1 - l)n + 7 2 + a/4(1 - 7)"' + 7 2 ) 



(135) 



where the central charge c and the parameter 7 are related by Eqs. (123). This dou- 
bly multifractal spec trum possesses the requested properties, like sup a //2(a, a') = 
f(a), where f(a) is ( |l22| ) above. 

This double multifractality can be generalized to higher ones by considering star 
config urati ons made of m simple scaling paths all originating at the same vertex, 
as in (101), with the following poly-multifractal results. The m-order case will be 
given by 



r m (ni,n 2 , ...,n m ) = x rn (ni,n 2 , —,n m ) - 2 

x m (m.nj, ...,n m ) = 2V \& m + U' 1 (m) + U' 1 (n 2 ) + ... + U' 1 (n m ) 



(136) 



Here A m is the quantum gravity boundary scaling dimension of the m-star S m made 
of m ( simple) scaling paths. According to the star algebra we have: 

A m = mU- 1 fa) = yC/- 1 (£ a ) = - (137) 

whre X2 = X is the boundary scaling dimension of a scaling path, i.e., a 2-star, 
already considered in Eq. ( jU7| ), and such that U^ 1 (x) = 1 — 7. We therefore 
arrive at a total (boundary) quantum scaling dimension 

i=l 

such that 

x m (ni,n 2 ,...,n m ) = 2V(S m ) . (139) 
Using the shift identity (see (|115|)) 



U- I {n) = l + \V- I {n)- 



where we recall that 



V~ x (n) = V 4 (l-7)« + 7 2 

we also have 



8=1 

The multiple Legendre transform ( |103[ ) is performed as above for the case m = 2. 
We have 

ttl = = 2F ' (,5m) £ = V ' {Sm) (n<)] ' ' (W1) 

so we get 



v /4(l- 7 )^+7 2 ^(O 

or, equivalcntly 



(142) 



inverted into 



^ (n,) = ^ , (143) 



n t = V I — I . (1441 



One gets from Eqs. ( |140|) and (|143|) 



HlH" 



This allows the simple expression of /„ 



f m ({a,}) = 2 - V(S m ) +Y, a > V ( — ) ■ 

i=l \ ai / 



(146) 



Reordering the S m terms with use of (110), and recalling identity (145) for 8 V 
finally gives after some calculations the explicit formulae 



/m ({«, = l,..., m }) =2 + 



7 



1 



2(l- 7 ) 



1-7 



m > 



7 \ - 

4(1-7) ti' 



with 



(147) 



m 1 



a, 



V4(l- 7 )^ + 7 2 \ 2 



^ j4(l- 7 K+7 2 



3=1 



Substituting expressions ( 123 ) gives in terms of c 

25 -c 1 



fm ({ai=l,..., m }) = 



12 

1 - c 



(1-7) 



24 



(148) 



(149) 



The domain of definition of the poly-multifractal function / is independent of c and 
given by 



m 

i-jE^>o, 



(150) 



as verified by Eq. (148). The two-sided case (134) above is recovered for m = 2, 
while the self-avoiding walk case ( |78| ) is recovered for 7 = —1/2, c = 0. Notice 
that the case / m= i(ai) corresponds to the potential in the vicinity of the tip of a 
conformally invariant scaling path, and differs from the usual f{a) — swp a , f?{a, a')) 
spectrum, which describes the potential on one side of the scaling path. 

We can also substitute equivalent "electrostatic" angles 9i — ir/cti to variables 
cti. This gives a new distribution: 



fm ({0i: 



1}) = fm ({C*i= 



7 



2(l- 7 ) 4(1-7)^ 



9 m 

_7_ 



(1-7) 



-771 1 1 



^ m 

2tt ^ 



(151) 



The domain of definition of distribution f m is the image of domain (150) in 9- 
variables: 



6, < 2vr . 



(152) 



i=l 



The total electrostatic angle is thus less than 2ir, which simply accounts for the 
electrostatic screening of local wedges by fractal randomness, as expected. 

Th e m aximum of f m or f m is by construction obtained for rii = 0, Vi = 1, ...,m. 
Eq. ( |148[ ) gives the values of singularity exponents &i at the maximum of f m ; 



7r ml 1 



J = ^{ 1 --) ,Vi = l,...,m, (153) 
corresponding to a maximum value of f m of f m : 

sup/ m = f m ({di=i,..., m }) = / m [{8<=i,..., m }) = 2 - 2F(A m ) 

' , 7 2 1 



2(1-7) 8(1-7) 



to 2 . (154) 



The interpretation of the poly-multifractal spectrum can be understood as follows. 
The probability P({oti}) = P({6i]) to find a set of local singularity exponents {oti} 
or equivalent angles {9i} in the m sectors of an w-arm star is given by the ratio 

P m ({a t }) oc RM{oi})/ R ^Pf n (155) 

of the respective number of configurations to the total one. We therefore arrive at 
a probability, here written in terms of the equivalent electrostatic angles: 



^({0,}) , (156) 

-l 

1 ' 



1 2 27T 



fm ({6i}) - f m ({6i}) = -^y— {^rn^ 



7 



2 



4(1 - T ) 



For a large scaling star, the dominant set of singularity exponents {en}, or wedge 
angles {9i}, is thus given by the symmetric set of values Jl53| ). 



IX Analysis of Multifractal Dimensions and Spectra 

Let us collect the results for the one-sided functions r(n), D(ri), and f(a). Each 
conformally invariant random system is labelled by its central charge c, c < 1. The 
multifractal exponents r(n) and generalized dimensions D(n) of a simply connected 
CI boundary then read explicitly: 



/x 1/ 25 -c /24n+l-c \ /iro , 
r(n)=-(n-l) + — 1 ' ( 158 ) 



r(n) 1 / /24n+ 1 -c \ , . 



1 - c 



24 

with a Legendre transform : 



/(«) 



25 



48 



2a 



24 



-a, 



(161) 



a G (J,+oo) 



It is interesting to note that the general multifractal function (161) possesses the 
invariance property (p3), since it also reads 



f(a)-a 



25 



24 



1 t 1 

1 - - 2a - 1 + 

2 V 2a -1 



(162) 



Notice that the generalized dimensions D{n) satisfy, for any c, r'Jn = 1) = D(n 
1) = 1, or equivalently f{a = 1) = 1, i.e ., Makar ov's theorem ]3q ], valid for any 



simply connected boundary curve. From ( 15£ ,16C ) we also remark a fundamental 
relation, independent of c: 



3-2D(n) = l/a = e/ir. 



(163) 



We also have the superuniversal bounds: Vc,Vn, h = D(oo) < D(n) < D(n*) = |, 
hence < 9 < 2ir. We arrive at the geometrical multifractal distribution of wedges 
9 along the boundary: 



/»-'(!)-? 



25 



12 9(2ir 



(164) 



Remarkably enough, the second term also describes the contribution by a we dge t o 
the density of electromagnetic modes in a cavity ]66]| . The simple shift in (164). 
25^25 — c, from the c = case to general values of c can then be related to results 
of conformal invariance in a wedge p7| . The partition function for the two sides of 
a wedge of angle 9 and size R, in a CFT of central charge c, indeed scales as [p8[ 



Z(9)RiR 12 0(2^-0) 



(165) 



Thus, one can view the c dependance of result (164) as a renormalization of the 
number of sites with wedge angle 9 along a self-avoiding scaling curve by a partition 
factor [i?(#)] -1 , representing the presence of a c-CFT along such wedges. The 
maximum of f(a) corresponds to n = 0, and gives the dimension Dep of the 
support of the measure, i.e., the accessible or external perimeter as: 



D E p = sup a /(a) = f(a(n = 0)) 

w 2(1-7) 2 24 
This corresponds to a typical exponent 



VT~c (V25 - c - vT^c) 



(166) 
(167) 



d = a(0) = 1 



1 



1 



VT~^ (%/25^c - VT~c) 



1 \12 

This corresponds to a typical wedge angle 



(3-2DEP)- 1 . (168) 



9 = ixja = tt(3 - 2D EP ) 



(169) 



In analogy to probability ( |155| ) for multiple angles, the probability P(a) to find a 
singularity exponent a or, equivalently, P{9) to find an equivalent opening angle 9 
along the frontier is 

P(a) = P(9) oc RfW-fW . (170) 
Using the values found above, this probability can be recast as (see also f3gf| ) 



P(a) = P{6) oc exp 



24 



In R Vl^c - 



(171) 



where 



1 7T 1 

The multifractal functions / (a) — a = f(9) — § are invariant when taken at primed 
variables such that 



2tt = 6 



(172) 



which corresponds to the complementary domain of the wedge 9. This condition 
reads also D(n) + D{n') = 2. This basic symmetry, first observed and studied in 
|60[| for the c = result of |E9| , is valid for any conformally invariant boundary. 
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FIG. 5. Universal multifractal exponents r(n). The curves are indexed by the central 
charge c, and correspond respectively to the same colors as in Fig. 7 below: (black: 2D 
spanning trees (c = —2); green: self-avoiding or random walks, and percolation (c = 0); 
blue: Ising clusters or Q = 2 Potts clusters (c = |); red: N = 2 loops, or Q = 4 Potts 
clusters (c = 1). The curves are almost indistinguishable at the scale shown. 
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FIG. 6. Universal generalized dimensions D(n). The curves are indexed by the same 
colors as in Fig. 7 below but are almost indistinguishable at the scale shown. 



The multifractal exponents r(n) (figure 5) or generalized dimensions D(n) (figure 
6) appear as quite similar for various values of c, and a numerical simulation would 
hardly distinguish the different universality classes, while the f(a) functions, as we 
shall see, do distinguish these classes, especially for nega tive n, i.e. large a. In Fig. 
7 are displayed the multifractal functions /, Eq. (161), corresponding to various 
values of —2 < c < 1, or, equivalently, to a number of components N € [0,2], and 
Q G [0, 4] in the O(N) or Potts models (see below). 
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FIG. 7. Universal harmonic multifractal spectra f(a). The curves are indexed by the 
central charge c, and correspond respectively to: 2D spanning trees (c = —2); self-avoiding 
or random walks, and percolation (c = 0); Ising clusters or Q = 2 Potts clusters (c = i); 
N = 2 loops, or Q — 4 Potts clusters (c = 1). The maximal dimensions are those of 
the accessible frontiers. The left branches of the various f{a) curves are largely indistin- 
guishable, while their right branches split for large a, corresponding to negative values of 
n. 



The singularity at a = |, or 9 — 2n 7 in the multifractal functions /, or /, corre- 
sponds to boundary points with a needle local geometry, and Beurling's theorem j58| 
indeed insures the Holder exponents a to be bounded below by \ . This corresponds 
to large values of n, where, asymptotically, for any universality class, 



Vc, lim D(n) = — . 

n — >oo 2 



(173) 



The right branch of / (a) has a linear asymptote 



lim f(a) I a = n* = -(1 - c)/24. (174) 



The limit multifractal spectrum is obtained for c = 1, which exhibits an exact 
example of a left-sided MF spectrum, with an asymptote / (a — > oo, c = 1) — > | 
(Fig. 7). It corresponds to singular boundaries where / (9 —> 0, c = 1) = 1 = -Dep, 
i.e., where the external perimeter is everywhere dominated by " fjords" , with typical 
angle = 0. It is tempting to call it the "Ultimate Norway" . 

The a — > oo behavior corresponds to moments of lowest order n — > n*, where 
D(n) reaches its maximal value: Vc, D(n*) — |, common to all simply connected, 
conformally invariant, boundaries. This describes almost inaccessible sites: define 
Af (H) as the number of boundary sites having a given probability H to be hit by a 
RW starting at infinity; the MF formalism yields, for H — ► 0, a power law behavior 



with an exponent 



Af(H)\ H ^ Q ^H-( 1+n "> (175) 



Some particular cases are worth considering. An Ising cluster possesses a multi- 
fractal spectrum with respect to the harmonic measure (c = i): 

1 , 7 



r(n) =_(n-l) + — (V48JT + T-7), (177) 



49 / 1 



/H=i( v 3 -^lJ-^ ^(i+oo), (178) 
with a dimension of the accessible perimeter 

D EP = sup a /(a, C =i) = H (179) 

The Q = 4 Potts model provides an interesting example of a left-handed multifractal 
spectrum (c = 1) 

r(n) =i( n -l) + V»-l, (180) 

/(«) =^( 3 -^3y)' aG(|,+oo), (181) 

with accessible sites forming a set of Hausdorff dimenson 

r> E p = sup Q /(a,c=l) = ^, (182) 

which is also the maximal value common to all multifractal generalized dimensions 
D(n) — —zrT(n). Notice that the external perimeter which bears the electrostatic 
charge is a simple curve, i.e. a curve without double points, a self-avoiding or 
simple path. We therefore arrive at the striking conclusion that in the plane, a con- 
formally invariant scaling curve which is self-avoiding has a Hausdorff dimension at 
most equal to Z?ep = 3/2 H]. The corresponding Q = 4 Potts frontier, while still 
possessing a set of double points of dimension 0, actually develops a logarithmically 



growing number of double points (69^ 



X Geometric Duality in O(N) and Potts Cluster Frontiers 

The 0(N) model partition function is that of a gas Q of self- and mutually- 
avoiding loops on a given lattice, e.g., the hexagonal lattice |is|| : Z (n) = 
J2g K^ B N^ p , with K and N two fugacities, associated respectively with the total 
number of occupied bonds Mb, and with the total number Mp of loops, i.e., poly- 
gons drawn on the lattice. For N G [—2,2], this model possesses a critical point 
(CP), K c , while the whole "low-temperature" (low-T) phase, i.e., K c < K, has crit- 
ical universal properties, where the loops are denser that those at the critical point 




The partition function of the Q-state Potts model on, e.g., the square lattice, with 
a second order critical point for Q G [0, 4], has a Fortuin-Kasteleyn representation at 
the CP: Zpotts = J2u{c) Q^ p , where the configurations U(C) are those of reunions 
of clusters on the square lattice, with a total number Mp of polygons encircling all 
clusters, and filling the medial square lattice of the original lattice [|l8|,[l7). Thus 
the critical Potts model becomes a dense loop model, with a loop fugacity N — Q 2 , 
while one can show that its tricritical point with site dilution corresponds to the 
0(N) CP @. The OiNf and Potts models thus possess the same "Coulomb gas" 
representations |l8|Jl7| , |70| | : 

N = = -2cos7T 5 , (183) 

with g G [1, |] for the 0(N) CP, and g G [±, 1] for the low-T 0(N), or critical Potts, 
models; the coupling constant g of the Coulomb gas yields also the central charge: 

c = 1 - 6(1 - gf/g. (184) 

Notice that from the expression ( |112| ) of c in terms of 7 < one arrives at the 
simple relation: 

1 = 1-9, 9 >1; 7=1-1/5, 9 <1- (185) 

The above representation for N = \[Q E [0, 2] gives a range of values —2 < c < 1; 
our results also apply for c S (—00, —2], corresponding, e.g., to the O (N S [—2, 0]) 
branch, with a low-T phase for g G [0, |], and the CP for g G [§, 2 ]. 

The fractal dimension -Dep of the accessible perimeter, Eq. (|l67| ), is, once rewrit- 
ten in terms of g, and like c(g) — c(g _1 ), a symmetric function of g 

D EP = 1 + \g-^{l - 9~ X ) + ^(1 - 9), (186) 

where $ is the Heaviside distribution, thus given by two different analytic expres- 
sions on either side of the separatrix g = 1. The dimension of the hull's frontier, 
i.e., the complete set of outer boundary sites of a cluster, has been determined for 
O(N) and Potts clusters (6|], and reads 

D n = l + \g-\ (187) 

for the entire range of the coupling constant g G [5, 2]. Comparing to Eq. ( |186| ), 
we therefore see that the accessible perimeter and hull dimensions coincide for 
g > 1, i.e., at the 0{N) CP (or for tricritical Potts clusters), whereas they differ, 
namely Z?ep < Dh, for g < 1, i.e., in the O(N) low-T phase, or for critical Potts 



clusters. This is the generalization to any Potts model of the effect originally found 
in percolation |64| . This can be directly understood in terms of the singly connecting 
sites (or bonds) where fjords close in the scaling limit. Their dimension is given by 

n 

^sc = i + \g~ x - \g- (188) 

Thus, for critical O(N) loops, g s (1,2] and Z?sc < 0, so there exist no closing 
fjords, which explains the identity: 

-Dep - Da; (189) 

whereas D$c > 0, g € [h, 1) for critical Potts clusters, or in the O(N) low-T phase, 
where pinching points of positive dimension app ear in the scaling limit, so that 



Dep < £>h (Table 1). we then find from Eq. (|186|), with g < 1 



(Dep-1) (£>h-1) = ~- (190) 

The symmetry point Dep = -Dh = § corresponds to g = 1, iV = 2, or Q = 4, where, 
as expected, the dimension -Dsc = of the pinching points vanishes. 

For percolation, described either by Q — 1, or by the low-T 0(N — 1) model, 
with g = |, we recover the result Dep = §, recently derived in p8[ . For the Ising 
model, described either by Q — 2,g — |, or by the 0(N — 1) CP, g' = g^ 1 = |, 
we observe that the EP dimension Dep = ^ coincides, as expected, with that of 
the critical 0(N = 1) loops, which in fact appear as EP's. This is a particular case 
of the geometric duality relation between the external perimeters of critical Potts 
models and the loops of O(N) models at their critical points: 

^ep (Q(g)) - D H (O (N(g'))) , for g' = g~ x , g<l. (191) 

In terms of this duality, the central charge takes the simple expression: 

c=(3-2 9 )(3-2 9 '). (192) 
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2 


3 


4 


c 


-2 





1/2 


4/5 


1 


Dep 


5/4 


4/3 


11/8 


17/12 


3/2 


Dh 


2 


7/4 


5/3 


8/5 


3/2 


Dsc 


5/4 


3/4 


13/24 


7/20 






TABLE I. Dimensions for the critical Q-state Potts model; Q — 0, 1, 2 correspond 
respectively to spanning trees, percolation and Ising clusters. 



XI Relation to the SLE K process 



An introduction to the stochastic Lowner evolution process (SLE K ) can be found 
in [fn| . This process essentially describes the cluster hulls we have introduced above. 
They can be simple or self-intersecting paths. The SLE K is parameterized by k, 
which describes the rate of an auxiliary Brownian motion, which is the source to the 
process. When k G [0, 4], the random curve is simple, while for k G (4, 8), the curve 
is a self-intersecting path. For k > 8 the path is space filling. The correspondence 
to our parameters, the central charge c, the string susceptibility exponent 7, or the 
Coulomb gas constant g, is as follows. 

In the original work by Schramm fl39| , the variance of the Gaussian winding angle 
of a SLE K of size R was calculated, and found to be V Klni?. In J73| we found, for 
instance for the O(N) model, the variance y/ (4/g)lni?, from which we immediately 
infer the identity 



4 

K = — 

9 



(193) 



The low-temperature branch g G [i, 1) of the O(N) model, for N G [0, 2), indeed 
corresponds to k G (4,8] and describes non simple curves, while N G [—2,0], 5 G 
[0, |] corresponds to k > 8. The critical point branch g G [1, |],7V G [0,2] gives 
k G [§, 4], while g G [|, 2], N G [-2, 0] gives k G [2, §]. The range n G [0, 2) probably 
corresponds to higher multicritical points with g > 2. Owing to Eq. ( |185| ) for 7, we 
have 



4 , _ K 

1 , «<4; 7 = 1 , k>4. 

k 4 



(194) 



The central charge (112) or (184) reads accordingly 



c = 1-24(^-1) Ik 



(195) 



an expression which of course is symmetric under the duality k/4 — » 4/k = k' , or 

kk' = 16 , 



reflecting the symmetry under gg' — 1 p4[ . The self-dual form of the central charge 
reads accordingly: 



C =i(6-«)(6-K')- 



(196) 



From Eqs. (187) and (186) we respectively find 



(197) 



Dbp = 1 + -i?(k-4) + £i?(4-k) , (198) 
k 8 

in agreement with later derivations in probability theory J74| . For k < 4, we have 
Dep(k) =_Dh(k). For k > 4, the self-intersecting scaling paths obey the duality 
equation (|190|) derived above, recast here in the context of the SLE K process: 



[£>ep(k) - 1] [D H ( K ) - 1] = -, k>4, 



(199) 



where now 



Dep(^) = D h {k' 



16/rc) k' < 4 . 



Thus we predict that the external perimeter of a self-intersecting SLE K >4 process 
is by duality the simple path of the SLE(i6/ K \ =K /<4 process. 

The symmetric point k — 4 corresponds to the 0(N = 2) model, or Q = 4 Potts 
model, with c = 1. The value k = 8/3, c = corresponds to a self-avoiding walk, 
which thus appears |?o|,[38[ as the external frontier of a k — 6 process, namely that 



Work remains to be done to elucidate the relation between the SLE construction 
in probability theory and the Coulomb gas and conformal invariance approaches, 
as well as the quantum gravity method described here. 
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